Angular Structure of Lacunarity, and The Renormalisation Group 
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We formulate the angular structure of Lacunarity in fractals, in terms of a symmetry reduction of 
the three point correlation function. This provides a rich probe of universality, and first measure- 
ments yield new evidence in support of the equivalence between self-avoiding walks and percolation 
perimeters in two dimensions. 

We argue the Lacunarity reveals much of the Renormalisation Group in real space. This is supported 
by exact calculations for Random Walks and measured data for percolation clusters and SAW's. 
Relationships follow between exponents governing inwards and outwards propagating perturbations, 
and we also find a very general test for the contribution of long range interactions. 
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In this letter we bring together two outstanding issues 
in the theory of fractals, which we believe will also have 
bearing on critical phenomena more generally. The first 
is 'Lacunarity' J^-Q] , which as originally introduced could 
entail all the discriminations of fractal structure (beyond 
dimension) which are evidently visible yet hard to codify 
in a simple generic way. The second is the Renormalisa- 
tion Group [0 which has been a predominant influence 
for three decades in the theory of critical scaling phenom- 
ena, fractals included, yet has never to our knowledge 
been regarded as a directly measurable object in its own 
right. Here we will present evidence and arguments that 
a particular version of Lacunarity measurement is effec- 
tively a measurement of the (Linearised) Renormalisation 
Group. 

The 'Lacunarity Function' was introduced || as a scal- 
ing reduction of the three point correlation function for 
a statistical fractal, 



not multifractal |6|,[7| any simple ambiguity in the mea- 
sure is expected to cancel from the Lacunarity Function 
C. 



I. THE ANGULAR STRUCTURE OF 
LACUNARITY 

We introduce here the decomposition of the Lacunar- 
ity Function over angular harmonics, for example in two 
dimensions 



C(R,0) =C (R) + 2 V C m (R) cosmO. 



(2) 



C(R,9) = 



C3(fi,r 2 ,r 3 ) 
C 2 (ri,r 2 )C 2 (fi,r 3 ) 



(1) 



Our physical interpretation of the individual components 
£m(f'/ r ) is that they measure the correlation between 
different lengthscales r and r' of relative fluctuations, of 
a particular angular harmonic, in the local two point cor- 
relation about each mass point. They can be more pre- 
cisely calculated as ratios 



where C 3 {fi,r 2 ,r 3 ) is the conditional probability for 
points r*2 and r 3 to be occupied given that fi is occu- 
pied, C%(f\,r2) the conditional probability for point r 2 
to be occupied given that fi is occupied, and the Lacu- 
narity Function C depends only on the angle 9 between 
the two vectors r 2 — fi and r 3 — r\ and their length ratio 
R = riz/r 12 . 

The Lacunarity function defined above is (by construc- 
tion) independent of absolute lengthscale and is a pure 
critical-point scaling object. It is crucial that we are dis- 
cussing measurements averaged over an ensemble of frac- 
tal objects, or at least self-averaging over the interior of 
one object much larger than than both r\ 2 and r\ 3 . The 
dependence of C on only a ratio of lengths follows from 
the assumption of continuous scale invariance, whilst the 
dependence on only a single internal angle follows in the 
more restrictive assumption that the ensemble of fractals 
in question has no favoured directions. The whole anal- 
ysis presupposes the definition of some suitable measure 
for the fractal structure in question, but provided this is 



C m (r'/r) = M m (r',r) 



< 5' m (r / ;ri)5 m (r;ri)* ; 

< S {r';n) >< S {r,fi) > 



(3) 



where S' m (r; fi) is the weighted count of occupied points 
within a shell of radius r about mass point r*i , each point 
weighted by the angular harmonic of its direction relative 
to fl (i.e. e lm6 in two dimensions), and the averages are 
with respect to all occupied mass points fi. jlC[ 

The width of the shells used to compute the averages 
in (S) drops out of the result provided it is small enough, 
but with finite datasets some compromise must be made 
to obtain reasonable signal to noise ratio. We typically 
used shells of 7% width by radius. In the case of diagonal 
elements (only), that is r — r' equivalent to R = 1, shot 
noise makes a systematically positive contribution to the 
reading for which we have corrected our data below. 

Figure 1 shows contour plots of the shell mass corre- 
lations M m (r',r) as a function of ln{r) and ln(r') for a 
large site percolation cluster grown at criticality. The 



1 



structure along the diagonal is direct evidence for the re- 
duction M m (r',r) = C m {R — r'/r) and on the basis of 
such plots we have averaged data parallel to the diagonal, 
within suitable looking windows, to obtain the curves for 
C m {R) shown in figure 2. 

Figure 3 shows a direct comparison between the corre- 
sponding measurements for the outer perimeter set of 
such a percolation cluster with those for a large self- 
avoiding walk. The broad agreement (with no adjustable 
parameters in any of the cureves shown) reinforces the 
claim that these two objects are in the same Universality 
Class, as very strongly suggested by Duplantier's caclu- 
lation H °f their identical f(a) spectra for the exterior 
harmonic measure. 

II. CONNECTION TO THE REAL SPACE 
RENORMALISATION GROUP 

The physical idea behind the renormalisation group 
is that under critical conditions the structure of a sys- 
tem is propagated through a continuous cascade of length 
scales. This is usually viewed in terms of the propaga- 
tion of certain (judiciously chosen) renormalised param- 
eters which remain invariant at the fixed point. The lin- 
earisation of the renormalisation group about the fixed 
point then corresponds to the transmission of perturba- 
tions through the lengthscale cascade. 

Our lacunarity function measures the correlation of 
structure fluctuations between different lengthscales and 
it is natural to suppose that this correlation arises due 
to propagation. If so, then our measurements of correla- 
tion directly probe the linearised renormalisation group. 
More strictly, we need the covariances rather than the 
correlations to probe propagation, to which end it is use- 
ful to distinguish from Lacunarity notation (and rescale) 
by writing 

cn m {R) = C r m{ * ) ~l ma , (4) 

where we propose to take the liberty of refering to C1Z 
as the Local Renormalisation Group Function. 

Our proposal is that if the local two point correlation 
function structure at radius r about point r\ is perturbed 
C2(r)p 1 — > C 2 (r)(l + e m (r;ri)H m ), where H m is a (nor- 
malised) angular harmonic, then the corresponding re- 
sponse at larger lengthscale r' should follow 

e m (r';n) = ClZ m (r' /r)e m (r; n). (5) 

It is important to note that the Local Renormalisation 
Group Function defined above is a 'shattered' response 
function, in that it gives the response due to changing the 
two point correlation at one locality. If we make a global 
relative perturbation e m (r) of the local two point cor- 
relation at lengthscale r around every mass point, then 
the growth of corresponding response defines a "Global 
Renormalisation Group Function" Q1Z through 



em{r') = gn m ir' /r)e ro (r). (6) 

From the number of (occupied) localities of radius r con- 
tributing to perturbation in a larger vicinity of radius r' , 
we expect the components of QTZ to scale as 

gn m {R) ~ R D az m iR) (7) 

for a simple mass fractal of fractal dimension D. Note 
that only even m are of physical consequence here, as the 
(global) two point correlation function is by construction 
an even function restricted to even perturbations. 

An important consequence of equation (ffl) is that 
global perturbations are relevant at large lengthscales if 
the corresponding Local RG function falls off more slowly 
than R~ D for large R. 

The case of simple Random Walks in d > 2 illustrates 
and tests these ideas with exact results. The general form 
of the Lacunarity function £(i?, 0) was first given in rcf 
H , but the angular harmonic decomposition of C1Z gives 
the much simpler form: 

OloiR) = R 2 - d ; OZmiR) = 1 + ^ - R~ m , m > 0. 

(8) 

This gives giZ (R) ~ R A - d , QK 2 (R) ~ R° and 
QTZ^iR) ~ R~ 2 for large R, which all have physical in- 
terpretation. 

A negative m = perturbation at lengthscale r means 
depressing the incidence of other points at radius r from 
any given one, corresponding in polymer language to 
turning on some excluded volume. That the resulting 
perturbation at larger lengthscales r' grows as (r'/r) 4_d 
corresponds to the well known Fixman perturbation ex- 
pansion of Coil Swelling ||]. The analogous positive 
perturbation leads to coil shrinkage (and ultimately col- 
lapse) . 

An m = 2 perturbation corresponds to the local struc- 
ture being stretched in some directions and shrunk in 
others, corresponding for a random walk to anisotropy in 
the individual walk steps. The scale independence of the 
resulting response follows from the equivalence of such 
random walks to affine deformations of isotropic ones, 
which we discuss as a much more general result below. 

An example of an m = 4 perturbation would be the 
intrinsic bias of random walks made on a square (more 
generally hypercubic) lattice, and the decay of the per- 
turbation reflects the known irrelevance of this in the 
large lengthscale limit. Analogous interpretations apply 
to higher orders of lattice bias, for example m = 6 and 
the triangular lattice in two dimensions. 

An important issue is whether the propagation of per- 
turbation is necessarily confined to running upwards in 
lengthscale. In systems governed by an equilibrium en- 
semble we can equally expect reverse (inwards) propaga- 
tion governed by the form of C1Z. Thus we are led natu- 
rally to predict that for every global perturbation propat- 
ing outwards from r to r' as (r'/r) _CT there is propagation 
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inwards of the same symmetry of perturbation from r' to 
r as (r' /r)~ T where r = a + D. 

The corresponding inwards perturbations of a random 
walk can be identified. The zeroth harmonic corresponds 
to simple confinement (and checks out), whilst higher 
harmonics can then be interpreted in terms of the influ- 
ence of a confining shape. The m — 2 case, ellipticity, is 
particularly important because this corresponds to poly- 
mer elasticity. 

Some of these interpretations can also be checked for 
less trivial fractals. For the Self Avoiding Walk (as for 
the random walk), the non-trivial part of the m = 
lacunarity is associated with the walk going from length- 
scale r out to lengthscale r' and coming back again. This 
is then governed by the probability of an SAW to close 
a loop, which is known to vary as N^ 1 oc r' where 
N is the number of steps in the walk (to reach radius 
r'). This leads to QlZo(R) ~ R° for large R, corre- 
sponding (correctly) to a change in the excluded volume 
causing the SAW to dilate uniformly on larger length- 
scales. The inset in figure 3 shows our measured data 
for Lo(R) — 1 oc CTZq(R) out to large R and is consistent 
with R~ D = i?~ 4 / 3 , albeit the data is noisy. The m — 2 
case for the SAW can be interpreted in terms of the scal- 
ing of elasticity for swollen polymers and leads to also to 
£R,2{R) oc R~ D at large R; this checks quite well against 
the main data in figure 3. 

We believe a very general result governs the m = 2 
lacunarity due to its RG connection, which explains why 
in all the fractals considered in this letter (numerically 
for the percolation cluster of figure 1, and analytically 
for random and self- avoiding walks), 

gn 2 (R) ~ r°. (9) 

This corresponds to the outwards global m=2 perturba- 
tion being marginally relevant, and it arises because a 
in = 2 perturbation can always be accommodated by 
a global strain of the metric of space. This should of 
course fail when long range interactions are important, 
sensitive to the metric of distance, and therefore equa- 
tion (||) serves as a general test of whether long range 
interactions play a role in the underlying physics. 

III. CONCLUDING REMARKS 

We have presented arguments that the Lacunarity 
Function £, or at least its asymptotes, can be viewed as 
response functions determined by the Renormalisation 
Group. It is harder to firmly establish just how much 
of the Renormalisation Group is thereby revealed. We 
can certainly restrict any such claim to the Linearised 
Renormalisation Group, because the Lacunarity Func- 
tion can only probe fluctuation correlations about the 
critical point behaviour: in RG terms this means that 
it can only probe the vicinity of the corresponding fixed 
point. We have restricted our attention to simple mass 



fractals, directly characterised in terms of one (positive 
definite) scalar field, the mass density. We speculate that 
C incorporates all the linearised RG relating to how this 
field might be coupled locally and bilinearly to itself, 
whereas it cannot know anything about what happens 
when other fields are introduced. 

A rich range of issues remain to be explored. We al- 
ready mentioned that the odd angular harmonics in the 
Lacunarity Function lack a simple Response Theory in- 
terpretation because the two point correlation function 
must be even. One way to deal with this is to recognise 
more explicitly the way the Lacunarity Function depends 
only on the shape of the triangle of points it correlates, 
which leads to a parameterisation of C involving only 
even harmonics; all the odd harmonics are thereby re- 
lated to the even ones. However it is more interesting to 
consider structures from irreversible growth where cor- 
relating present with future growth is the natural cor- 
relation to investigate and odd harmonics are no longer 
forbidden in such time resolved correlations. Wc antici- 
pate that problems like the response of Diffusion Limited 
Aggregation to anisotropy should be analysed in this way. 

In principle our work opens the door to measuring 
the Renormalisation Group, via Lacunarity, directly from 
real experimental data (as opposed to just computer sim- 
ulations). We recognise, however, that this may not be 
easy. Our simulations span some three decades of length- 
scale (i.e. cluster radii of order 1000 units), and even 
then the match up to expected scaling laws is not always 
good. The best way forward would appear to be not to 
target the asymptotic exponents, but rather to compare 
absolute lacunarity values between experiment, simula- 
tion and theory at moderate R. On the theory side, our 
exact Random Walk results are extendable to Markov 
Chains, but no other exact results are known (other than 
exponents). 

Finally, whilst all the fractals discussed in this letter 
have close relation to equilibrium critical phenomena, our 
present definition of Lacunarity certainly does not extend 
trivially to critical phenomena in general. We believe the 
correct analogue of our three point correlation for a mass 
fractal is a four point correlation of the corresponding 
critical distributions, the extra point being at infinity. 
The natural generalisation of our lacunarity function to 
spin models would thus be in terms of a four point corre- 
lation function at the critical point, and we look forward 
to exploring this in future work. 
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FIG. 1. Contour plot of the cross correlation of local two 
point function between different lengthscales, as per equa- 
tion ([]) with m — 0, for a percolation cluster of 1.5 10 6 sites 
grown on a traingular lattice at the percolation threshold. 
The axes are labelled proportional to the logarithm of the two 
lengthscales correlated. Contour lines parallel to the diagonal 
correspond to this reducing to a scale invariant 'Lacunarity 
Function' depending only on the lenghtscale ratio. Scale in- 
variance breaks down around the edges of the plot where one 
lengthscale approaches either the lattice spacing or the sample 
size. 
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FIG. 2. The different angular harmonics of the Lacunarity 
Function for the two dimensional percolation cluster of figure 
1. The m = data correspond to the central part of the 
contour plot averaged parallel to the diagonal, and the other 
harmonics are defined analogously in equation (bl) 




2.0 



m=2 



m=0 











\\ "--a 




- \ 




\ ■ "••■a. 


° m=0 


A , 
" \ " 








\ 




X \ 




\ \ 




\ 








m=4 






m=2 


m=6 





-10.0 



0.0 1.0 2.0 3.0 4.0 5.0 

-12.0 

FIG. 3. Angular harmonics of the Lacunarity Function 
compared between the Outer Perimeter of a Percolation Clus- 
ter (as per figure 1) and a Self- Avoiding Walk (in two dimen- 
sions) of comparable size. The general agreement of these 
measurements, which contain no adjustable parameters, sup- 
ports both types of object being in the same unversality class. 
The inset shows the data for the zeroth harmonic extended 
to larger lengthscale ratios. 
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